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Abstract. We present an analysis of high precision V light curves (LC) for 18 local Type la 
Supernovae, SNe la, as obtained with the same telescope and setup at the Las Campanas Obser- 
vatory (LCO). This homogeneity provides an intrinsic accuracy a few hundreds of a magnitude 
both with respect to individual LCs and between different objects. Based on the Single De- 
generate Scenario, SD, we identify patterns which have been predicted by model calculations 
as signatures of the progenitor and accretion rate which change the explosion energy and the 
amount of electron capture, respectively. Using these templates as principle components and the 
overdetermined system of SN pairs, we reconstruct the properties of progenitors and progenitor 
systems. All LCO SNe la follow the brightness decline relation but 2001ay. After subtraction of 
the two components, the remaining scatter is reduced to ~ 0.01 to 0.03^. Type SNe la seem 
to originate from progenitors with Main Sequence masses, Mms > 3 Mq with the exception of 
two subluminous SNe la with Mms < 2Mq. The component analysis indicates a wide range of 
accretion rates in the progenitor systems closing the gap to accretion induced collapses (AIC). 
SN1991t-like objects show differences in dml5 but no tracers of our secondary parameters. This 
may point to a different origin such as DD-Scenario or the Pulsating Delayed Detonations. 
SN2001ay does not follow the decline relation. It can be understood in the framework of C-rich 
WDs, and this group may produce an anti-Phillips relation. We suggest that this may be a result 
of a common envelope phase and mixing during central He burning as in SN1987A. 



1. Introduction 

SNe la are thermonuclear explosions of WD ( |Hoy lefcFowler , 1 96Q| ) i.e. end-stages of 
stellar evolution of stars between 1 and 7 Mq. Most likely, they result from the explosion 
of a C/O-WD with a mass close to the Chandrasekhar limit (Mch)^ which accretes matter 
through Roche- lobe overflow in a single degenerate scenario (SD) ( [Whelan fc Iben,1973^ ), 



or merging of two degenerate WDs (DD) ( |Webbink,1984[ [iben fc Tutukov,1984 ). We re^ 



gard SDs as most likely for the majority of SNe la because e.g. of the homogeneity in 
LCs and spectra, though there is strong evidence of contributions of both to the SNe la 
population (see |Hoefiich,2QQ6] and references therein). One of the keys is the empiri- 



cal relation between maximum brightness and the rate of decline, drnlS, (Phillips, 1993). 



From theory, dml5 is well understood: LCs are powered by radioactive decay of ^^Ni 
( [Colgate fc McKee,1969D . More ^^Ni increases the luminosity and causes the envelopes to 
be hotter. Higher temperature means higher opacity and, thus, longer diffusion time scales 
and, thus, slower decline rates after maximum light ( Hoefiich et al., 1996, Nugent et al.,1997[ 



Maeda et ah, 2003} [KasenfcWoosley,2QQ9 ) . The existence of a cZml5-relation holds up for 
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Longer diffusion TS 



Change by expansion rate/geometrical dilution 



Figure 1. Principle components based on theoretical models. We show the difference in V 
brightness Am(t) (in magnitudes) as a function of time (in days) relative to a reference model 
with solar metallicity and a main sequence mass Mms of 7Mq and a central density of the 
exploding WD of pc of 2 x lO^^f cm~^ . The plots are for models with TM© and 5 x lO^^f cm~^ , 
respectively. The annotation on the graphs give the main reason for the differences. 



virtually all scenarios as long as there is an excess amount of stored energy to be released 
dHoefiich et aT^lQ QG ) . 

Within SDs, the favorite models are Mch explosions in which burning starts as de- 
flagration which, at some point, transitions to a detonation DDT ( Khokhlov,1991) . The 
^^Ni production depends mostly on the DDT, i.e. one parameter. The expanding SN- 
envelopes have similar velocity and density structures because same masses, most of the 
WD undergoes burning, and the nuclear binding energies of the burning products are 
nearly the same. A dispersion of 0.2. ..0.3^ is expected (see Fig. [1]). Its origin can be 
related to the progenitor and accretion from the donor star: Changes in the metallicity 
and the Main Sequence Mass, Mms^ will change the size of the explosion energies by 
about 20 % because their influence of the size of C-depleted core formed during stel- 
lar He-core burning (Hoefiich et al.,1998, Dommguez et al.,2001[ ). Increasing accretion 
rates will decrease the central WD density pc at the time of explosion and, consequently, 
an increased electron capture will produce more stable iron-group elements at the ex- 
pense of ^^Ni. High precision LCs have been obtained at LCO for 18 local SNe la 
( [Contreras et al.,201Q^ Folatelli et al.,2010) and the individual, theoretical signatures 



(0 



have been recovered ( [Hoefiich et al.,20TQ ) . Note that a similar pattern as by pc can be 



produced by off-center ^^Ni distributions but Am increases later at ^ 60 to 100 days. 



2. Diversity of Type la Supernovae 

Secondary parameters: Two SNe la may differ in both pc/accretion and the C- 
depleted core/ Mms- Therefore, we employed component analyses to study secondary 
parameters. We use V because this color is hardly effected by metallicity, asphericity 
effects, and k-corrections. The differences Am{t) in LC pairs are described by 

Arriij^obsit) = ^ ><k{ij)fk{t) ^ Resij{t) 

k=l,2 

with Xk{ij) being the coefficients for a pair of SN i and j , fk{t) the principle compo- 
nents, and Res{t) the residuals. In our sample, we include 18 SN from LCO, i.e. 153 
pairs. Only 36 Xij = gk{i) / dk^jY^ are independent where gk{i) is the eigenvalue of fk 
to be attributed a specific SNe la. By solving the overdetermined system for gi us- 
ing a Simplex Method ( [NeldefcMead, 1965| ) , we obtain most likely values for Xk{ij). In 
Fig. [21 the pairs SN2005al/SN2005am, SN2005el/sn2005ef, SN 2005ef/SN 2005na and 
SN 2005al/SN2005ef are given for A(ij). Overall, the residuals are consistent with zero 
but fits are not unique due to errors in brightness and time coverage (see Fig. [3]). Remap- 
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c) very similar or method does not work 




Figure 2. Best fits to the observations of individual pairs of SNe la. We give tfie weigfited 
components, their sums along the observations (with error bars) and the residuals (crosses) for 
cases which are dominated by the central density /accretion rate (upper left), in the progenitors 
(upper right), similar pc and M© (lower left), and a mixed case (lower right). Note that the 
residuals are small and, within the error bars, consistent with zero. 




Figure 3. Probability distributions of Xij for the pairs of SN2005al/SN2005am and 
SN2005el/sn2005ef based on MC solutions for the overdetermined system. Sparse time coverage 
produce large uncertainties in the eigenvalues. 



ing the individual gk{i) to Mms and pc, shows that (a) pc are evenly distributed from 
1 X 10^^ cm~^ to 7 X 10^ gcm~^, i.e. close densities leading to an accretion induced 
collapse (AIC) (b) SN la come from massive progenitors with M ^ 3Mq but the two 
subluminous SNe la indicate Mms between 1 ... 2 Mq. All pairs of SN1991t-like objects 
show Xk{ij) = in all components. Either, they are very similar despite differences in 
(iml5, or they lack a central region of high densities and similar C/0 in the progenitors 
as may be expected for mergers (DD) or pulsating delayed detonation models. For more 
details and a complete analysis of all 28 LCO SNe la, see Sadler et al. (in preparation). 
SN2001ay shows that nature is even more diverse. SN2001ay is slower than any SNe la 
known, combined with a fast rise of some 16 days (Krisciuanas et al. 2011). SN2001ay 
would be brighter by about 1^ based on dm 15 and the distance to the host galaxy. In 
fact, dmlb is slower than implied by the instantaneous energy input by radioactive decays 
Ej. We submit that, still, this SN can be understood within the physics underlying the 
dml5 relation, and in the framework of pulsating delayed detonation models originating 
from a Mch mass WD but with a core of some 80% C rather instead of the 15 to 20 
% usual for stellar central He burning. Higher C fraction means more nuclear energy by 
^^C{a^ 7)^^0 by ~ 40 % and faster expansion of the inner layers. Faster expansion means 
that a larger fraction of the energy by ^^Ni decay goes into expansion work rather than 
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SNSOOlay vs. PDD _llb 




Figure 4. B and V LCs of SN2001ay (Krisciunas et al.,2011) in comparison with theory. We 
give the instantaneous deposition by radioactive matter for FDD lib, a pulsating model with 
a carbon rich core, and a classical DD model (upper left), and the B and V (lower left). The 
DD model has been scaled by 0.037 dex to equal the gamma-ray input at maximum light. The 
comparison between SN2001ay and FDD lib in B (upper right) and V (lower right) as a function 
of time since maximum light in V observed. We assume a distance module m-M of 35.75"^, and 
a total reddening E(B-V) of 0.06"^. 

boosting the luminosity at maximum light, lower optical depth and shorter rise time. In 
our models, the the maximum brightness is smaller than the instant radioactive energy 
release ( [Arnett, 19 88), and the light curves approaches "from below". Our model 
agrees reasonably well with the observations (Fig.|4]). The reason for a high C abundance 
we can only speculate. During the early stages of central He burning, high C abundances 
are produced by ^He{2a^ lY'^C burning but, at the end, ^^C(a, 7)^^0 depletes unless 
strong mixing of He avoids this phase. A possible path may be a common envelope like 
in the progenitor of SN1987A. For more details, see Baron et al. (in preparation). 
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